Decoherence in adiabatic quantum evolution — application to Cooper pair pumping 
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One of the challenges of adiabatic control theory is the proper inclusion of the effects of dissipation. 
Here, we study the adiabatic dynamics of an open two-level quantum system deriving a generalized 
master equation to consistently account for the combined action of the driving and dissipation. 
We demonstrate that in the zero temperature limit the ground state dynamics is not affected by 
environment. As an example, we apply our theory to Cooper pair pumping which demonstrates the 
robustness of ground state adiabatic evolution. 



Accurate control of quantum systems has been one of 
the greatest challenges in physics for the last decades. 
Adiabatic temporal evolution [1] has attracted a lot of 
attention [2-5] in this respect since it provides robustness 
against timing errors and typically utilizes evolution in 
the ground state of the system. Such evolution has been 
argued to be robust against relaxation and environmental 
noise [6-8]. 

The combined effect of adiabatic evolution and dissi- 
pation were considered by many authors using various 
techniques and with different aims and assumptions, see, 
e.g., Refs. [6, 9-13]. We derive in this Letter a unique 
master equation that treats the combined effect of noise 
and adiabatic driving consistently and, thus, provides a 
pioneering tool for studying the effects of decoherence in 
quantum control protocols employing adiabaticity [3, 4]. 
We find that adiabatic evolution should not be treated in 
the secular approximation [14]. Furthermore, the master 
equation incorporates new terms ensuring relaxation into 
the correct time-dependent ground state. When these is- 
sues are properly addressed, the expectation values of 
physical observables in the adiabatically steered ground 
state are not influenced by zero-temperature dissipation. 
We apply our theory to adiabatic charge transport in 
superconducting circuits in the presence of noise. In 
spite of its long history [15-18], this problem has re- 
cently attracted revived theoretical [19-22] and experi- 
mental [23, 24] interest due to its fundamental relation 
to geometric [25] and topological [26] phases and to its 
potential applications in metrology [24, 27]. 

We consider an open quantum system subject to exter- 
nal time-dependent control fields. The total Hamiltonian 
of the system and its environment, H(t), is the sum of 
three terms, H(t) = Hs(t)+HE+V, where Hs(t) denotes 
the time-dependent system Hamiltonian, He is the bath 
Hamiltonian and V is the system-bath coupling. Assum- 



ing that the driving does not directly affect the coupling 
term between the system and the environment, we can 
write V = X ® Y, where A is a bath operator and Y 
is a system operator. In the case of weak system-noise 
coupling and slow driving, a convenient basis to describe 
the dynamics of the system is the instantaneous energy 
eigenstate basis, also called adiabatic basis, defined by 
H S (t) |Vn(*)> = E n {t) IV'n(t)). The states \if> n (t)) are as- 
sumed to be normalized and non-degenerate. We denote 
by D(i) the transformation from a given fixed basis to the 
adiabatic one. The evolution of the transformed density 
matrix is governed by the effective Hamiltonian 



ffW(t) = H s {t) + hw(t)+V(t)+H E 



(1) 



where H s (t) = (t)H s (t)D(t), V(t) = D^(t)VD(t) = 
X (g) Y(t), and w = —iD^D. 

We note that there are a few possible strategies of 
treating the dissipation. The usual one is to disregard 
w in the calculation of the dissipative rates [6]. Then, 
zero-temperature environment tends to relax the system 
to the ground state of Hs(t), while the rotation w tries 
to excite the system. The resulting state is different from 
both the adiabatic ground state (ground state of Hs) and 
from the ground state of Hs + Hw. The second strat- 
egy is to first perform a series of transformations to the 
super-adiabatic bases [13, 28] and, then treat the dissi- 
pation. The first step would be to diagonalize Hg + Hw 
with a unitary transformation D± and get a much smaller 
non-adiabatic correction w\ = —iD\t)\. Here, the dis- 
sipation (treated in Markov approximation) takes us to 
the ground state of Hs + hw. Although not exact, the 
second strategy allows one to treat the combined effect of 
noise and driving consistently. Here we adopt this strat- 
egy to calculate the lowest order correction to the adia- 
batic dissipative dynamics of a two-level system. As we 
will show, up to higher order corrections, this treatment 
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correctly accounts for the relaxation to the ground state 
of the superadiabatic hamiltonian Hs + hw. By using 
standard methods explained, e.g., in Rcf. [14], we arrive 



at the following master equation for the reduced system 
density matrix pi(t) in the interaction picture (for the 
derivation see Appendix): 



dpi{t) - i[p I (t),w I {t)]-^ E { j cW [,>,(/) ^-.f/f/T.f/Cn 



dt 



n? 

dt' / dt" 



o jo 



[p I (t)®p E ,[w I (t , ),V I (t'')]]Mt) 
I 



(2) 



where Tr E indicates trace over the environmental de- 
grees of freedom and p E is the stationary density op- 
erator of the environment. To obtain Eq. (2) we have 
to take consistently into account corrections up to the 
order wVV, resulting in a nonstandard commutator ex- 
pression. The interaction picture operators are defined 
as Oj(t) = e iHB ^ h Ul(t,0)d{t)U s (t,0)e- iHEt / h , where 
U s (t,0) = e- l ti H s {T)dT/h ig the systcm time-evolution 
operator. In Eq. (2), the first contribution on the r.h.s. 
is of order a — h/( AT p ) where A is the minimum gap in 



the spectrum of Hs and T p is the period on which the 
Hamiltonian is varied [29]. The second term is as in the 
standard Bloch-Rcdfield theory The third one is a cross- 
term of the drive and dissipation ensuring relaxation to 
the proper ground state [13]. 

We now focus on the case of a general two-state system, 
with the instantaneous eigenstates \g) (ground state) 
and |e) (excited state). In this case, returning to the 
Schrodinger picture, we can recast Eq. (2) into 



and 



Pgg = - 2 ^Mw* ge Pge) - i^ge + ?eg)Pgg + T eg + f ^c( Pge ) 

+^[(2f + - f o)( i - Pgg) - (2 f_ - f )p„] + 2 ^^\ r ge + r eg - r o) , (3) 



p ge = iw ge (2p gg - 1) + i{w ee - w gg )p ge + iuj p ge - i(T ge + T eg )^sm(p ge ) - T v p ge + (T + + r_)p. 



gg - 1 + 



+ (2r _ _ T) _ .^mA^ Tge)]pgg _ [!V (2r+ _ T ) + ; Hv) (reg _ Tge)] (1 _ pgg) 

OJq OJq U!q U!q 

+2 [<l^ e (p ge ) + 2i ^£il 3m(p se )] (f - f + - f _). (4) 



By Oki we denote the matrix elements (ra|0|n) of a 
general operator O, with m,n = e, g, except w mn = 

Y 2 

—i(m\h). We have defined the rates T ge = -^-S{—ujq) 

Y 2 

(excitation), T eg — -j^S(+u>o) (relaxation), T v = 

Y 2 

2-ff-S(0) (dephasing), and the less common transition 
terms f± = S (±u; ) , f = 2^^S(0), I± = 

^fS(±LJ ), and T = 2^5(0). Here, the matrix 
elements of Y obey Y gg (t) = -Y ee (t) and Y eg (t) = 
Y ge (t) [30]. The energy separation between the two 
states is hu n , which varies along the pumping trajec- 
tory. The power spectrum of the noise is defined through 
S( W ) = /~^X 7 (r)X 7 (0)>e*-dT. 



Throughout, we have used Markov approximation, i.e., 
we neglect the variation of pi(t) between t and t + t c , as- 
suming that the correlation time of the bath, r c , is much 
shorter than the typical relaxation time of the system, 
1/r. Furthermore, we made the approximation of adi- 
abatic rates (AR), i.e., in the calculation of the rates 
we neglect the slow variation of cjo, Y, and w, assum- 
ing the bath correlation time to be much shorter than 
the driving period t c <C T p . On the other hand, Eqs. (3) 
and (4) include all the non-secular terms traditionally ne- 
glected [14]. They introduce cross-dependence between 
p gg and p ge in the dissipative terms, and in our problem, 
omitting them would lead to unphysical results, such as 



violation of charge conservation. 

We are interested in the quasi-stationary limit that 
the system reaches when the evolution is adiabatic and 
it is initially in the ground state. We thus look for the 
solutions of p gg = and p ge = for a -C 1. Since 
Wmn = 0(a), in the absence of dissipation, we find 
that pgg ~ 1 + 0(a 2 ) and p ge ~ ~w ge /iUQ + 0(a 2 ) 
are the desired solutions. In the zero-temperature limit, 
S(— ojq) = 0, to the first order in a, Eqs. (3) and (4) yield, 
again, p gg = l + 0(a 2 ) and the following equation for the 
off-diagonal element up to order a: iojQ^l ge — T v fl ge — 
iT eg Qm(flg e ) — 0, with Q ge = p ge + w ge /ujo. The solu- 
tion of this equation is exactly the same as for the closed 
system; p ge = —w ge /ujQ. Therefore, the ground state evo- 
lution is not influenced by coupling to a zero-temperature 
Markovian environment in the adiabatic limit. Note, that 
including the imaginary part of the rates, e.g., the Lamb 
shift, does not change this result. 

The vanishing of the effects of dissipation is consis- 
tent with the following simple argument. In the zero 
temperature limit, and to first order in a, the effect of 
dissipation is to bring the system to the instantaneous 
ground state of the effective Hamiltonian, H\ = Hs+hw, 
which means that in the eigenbasis of Hi spanned by 
the eigenvectors, l^ri )j the density matrix has the form 
Pmn = {4>m\p ) = S mg S ng + 0(a 2 ) independent 
of the dissipative rates. Thus, within our approxima- 
tions, the ground state evolution is robust against zero- 
temperature environmental noise and the expectation 
value of any operator in the quasi-stationary evolution 
does not depend on the specific properties of the envi- 
ronment. If, instead, we neglect the non-secular terms, 
we obtain the same solution for p gg but the evolution of 
Pg e is influenced by the noise as p ge = —w ge / (u>q + iT/2), 
where T represents a combination of the dissipative rates. 
This leads to different expectation values of physical ob- 
servables that depend on p ge , and to the loss of robust- 
ness of the ground state dynamics. Therefore, in general, 
the non-secular terms cannot be neglected: they give a 
leading order contribution in Ta/A to the dynamics. 

To test our theory on a concrete example, we discuss 
a superconducting Cooper pair pump. It consists of an 
array of Josephson junctions coupled to two supercon- 
ducting leads, being subject to time-dependent external 
fields. As discussed by various authors (see, e.g., Ref. 
[19]), the transferred charge is the sum of a dynamic and 
a geometric contribution, Q = Q D + Q G . The first one 
corresponds to the average supercurrent and the second 
one to pumping. Assuming that only two levels are in- 
volved, the two contributions to the charge transferred 
through junction i in a pumping cycle can be written as 

d f Tp 

Q? = / (Pggh.gg + Peeh,ee)dt, (5) 
JO 

Qf = f P 2$te(p ge I ijeg )dt, (6) 
Jo 
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FIG. 1: An example of a quantum pump, the Cooper pair 
sluice, is shown on the left. A pumping cycle is sketched on 
the right. The time-dependent classical control parameters 
are the magnetic fluxes tuning the Josephson tunnel couplings 
Jl and Jr and the gate voltage controlling the offset charge 
n g of the island. They vary in time with period T p , whereas 
the phase difference across the device, tp, is stationary. 

where Ii is the current operator through junction i. Here 
we focus on the pumped charge, i.e., Qf = f^lfdt [31]. 
By substituting p ge = ~w ge /ujQ in Eq. (6) we arrive 
at the well-known formula for the adiabatically pumped 
current in a closed system [18], if = — -£-$le{w ge Ii teg ). 
As discussed above, this is also the limit of the adiabatic 
evolution in the presence of environmental noise. 

In particular, we consider the Cooper pair sluice [27] 
of Fig. 1. It consists of a single superconducting island, 
coupled to two superconducting leads via two SQUIDs, 
i.e., Josephson junctions whose critical currents can be 
tuned by magnetic fluxes. The electrostatic potential 
on the island can be controlled by a gate voltage, V g , 
and there is a constant superconducting phase difference, 
(p = (fL — ipn between the two leads. In the absence of 
noise, the Hamiltonian of the sluice can be written as 

H.s = E c (n-ng) 2 -J L cos((p L -6)-J R cos(6-ip R ). (7) 

Here 9 and n are the operators for the superconducting 
phase of the island and the number of excess Cooper pairs 
on it. The Josephson couplings to left and right lead are 
denoted as Jl and Jr, n g = C g V g /2e is the normalized 
gate charge, and Eq — 2e 2 /Cs is the charging energy 
of the sluice; C g is the gate capacitance and Ce the to- 
tal capacitance of the island. The current operators of 
the left and right junctions read II = ^Jz,sin(^i — 6) 
and Ir = ^Jflsin(# — ifa), respectively. For Eq 3> 
max{ Jl, Jr}, and n g ~ 1/2 only two charge states, |1) 
and |0), i.e., one or no extra Cooper pairs on the is- 
land, are relevant. Dissipation is then mostly due to gate 
voltage fluctuations. Other noise sources, not considered 
here, are determined by fluctuations of the fluxes in the 
SQUIDs or in <p [19]. In the two-level approximation, the 
coupling between sluice and charge-noise has the form 
V = —g<J z ® fiV g (t), where g — eC g /C^ is the coupling 
constant, <r z = |0)(0| — |1)(1| and SV g (t) is the gate volt- 
age fluctuation. In the absence of dissipation, for the cy- 
cle Of Fig. 1, with Ji g [J mi „, Jmax], n g € [u g mjn, Tl g roax ] 
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FIG. 2: Pumped charge of the sluice under gate charge 
noise, (a) Phase dependence of the pumped charge Q G in 
a fully symmetric pumping cycle with respect to Jl and Jr 
of Fig. 1 as a function of the phase bias if. The dashed 
line shows the analytic result of Eq. (8) for adiabatic pump- 
ing. The solid lines are from the numerical calculations 
based on Eqs. (3) and (4) for / = T~ l = 100 MHz, with 
Cg/Cs = 0.015,0.0175,0.02,0.025, and 0.3 from bottom to 
top. (b) Coupling dependence of the pumped charge at 
ip = 7r/2. The dashed line shows the analytic result as in (a). 
The solid lines are for / = 10, 100, 150, 200 and 300 MHz from 
top to bottom. The other parameters are J max /-Ec = 0.1, 

Jmin/Jmax = 0.03, Tig max = 0.8, U g min = 0.2, E C / k B = 1 K 

(Ec/2irh = 21 GHz), R — 300 kQ, environment temperature 
T = 0, S(ljo) = 2hcj R, 5(-w ) = 0, and 5(0) = 2k B T R 
with Tb = 0.1 K. 



and for J max <C Eq, one obtains the pumped charge in 
the adiabatic limit according to Eq. (6) as 



2e(l 



J 11 



■ COS ip) 



(8) 



for both junctions [27]. Thus the transported charge 
depends on if, the average being one Cooper pair per 
cycle. In the presence of dissipation, Eqs. (3) and 
(4) were integrated numerically to obtain the temporal 
evolution of the density matrix along a pumping tra- 
jectory of Fig. 1. Figure 2 shows that, upon increas- 
ing the system-environment coupling at finite frequencies 
/ = Tp 1 , the pumped charge approaches the analytic re- 
sult of Eq. (8) for adiabatic pumping at all values of if, 
see Fig. 2(a). Figure 2(b) shows the coupling dependence 
of the pumped charge at various frequencies for if = n/2. 
On lowering the frequency, all the data collapse towards 
the horizontal dashed line which is again the result of 
Eq. (8). For / = 10 MHz, the numerical and analytic 
results are indistinguishable on this scale. Thus cou- 
pling to zero-temperature Markovian environment seems 
to be useful for adiabatic ground state pumping. We 
note, however, that Eqs. (3-4) are strictly valid only for 
adiabatic evolution and weak coupling. 



In conclusion, we derived a master equation for an 
adiabatically driven two-level system including the com- 
bined effect of drive and relaxation. We found it impor- 
tant to account for the time-dependence of the Hamilto- 
nian of the system in determining the dissipative rates 
and to include the non-secular terms. As an example, we 
analyzed adiabatic Cooper pair pumping in the ground 
state and demonstrated that the pumped charge is not 
influenced by zero-temperature environment. Numerical 
solution of the master equation suggests that dissipation 
can resume adiabatic pumping at finite frequencies. 
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Derivation of the master equation 

The master equation (2) is obtained with a develop- 
ment taking into account all the terms up to the order 
wVV. As in standard derivation of the master equation 
[1], we assume that the density matrix of the environ- 
ment is stationary and that the average of V over the 
environment degrees of freedom vanishes. 

Denoting the total density matrix of the system and 
environment as p tot (t) and employing the transforma- 
tion to the adiabatic basis as pj ot = D^p tot D, the von 
Neumann equation in the interaction picture reads 



p t I ot (t) = -[pf\t),hw I (t) + V I (t)]- 



(9) 



Notice that, in the weak coupling and in the adiabatic 



J 



limit, wi(t) and V/(t) are perturbative contributions of 
different order. Tracing over the degree of freedom of the 
environment, Eq. (9) becomes 



h (t) = i [pj (t) , wj (i)] + j:Tr E { IpT (t),Vi(t)]}, (10) 

where p\ ot = Tr E {pT}- 

Together with Eq. (9), we employ the identity 



( T ) = pf\ Tl ) + / dr'pr{r'). (11) 

Using iteratively Eqs. (9) and (11), we can obtain a per- 
turbation expansion of Eq. (10). 

Substituting Eq. (11) in the last term in Eq. (10) we 
have 

fc{t) = i[pi (t),wi(t)] + 

Since the average of V over the environment degrees 
of freedom vanishes, Tr^j [pj ot (0), V/(f)]} = and using 
Eq. (9), we obtain 



h{t) = i[p!(t), wj(t)} - ^Tt e { J dt'[[pf\t'), h Wl (t% Vj{t)] } - ^e{J^ dt'[[pr(t'), Mt'^Mt)] }. (13) 



The term /5j ot (i') can be transformed using Eq. (11); in the second and in the third term on the 
in particular, we substitute /5j ot (t') = /5/ ot (0) + right, respectively. Consistently, we have that 
/„*' df'pf^t") and pf^t') = ^(t) - f t , df'pf^t") r ^E{[[p t I ot (0),w I (t%V I (t)]} = and hence 



Pi{t) = i[p/(*W(*)] - ^Tt e { dt'[[pf\t)Mt%Vi{t)]} - ^T:e{ f & j* dt"[[ft°\t'%Hw I (t%V I (t)]} 
+ f Q dt ' f dt " Wit"), Vjit% Vjit)] }, (14) 



r 



Using Eq. (9) we eliminate pj ot (i") from the above obtain 
equation and, keeping the terms up to order wVV , we 
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Mt) = t [p I (t),w I (t)]-^Tr E { j\t'[[pf\t)Mt%Vi{t)}}-^TE{J^dt' dt'\[\? I °\t")Mt'%Mt%Vi{t)]} 



+ ^ TlE i f dt ' f dt"[[W\t"), wjimVj^)] , Vj(t)]}, 



(15) 



The third and fourth terms on the right are both of 
order wVV, namely, the highest order in our expansion. 
The last step in our derivation is to use Eq. (11) to sub- 
stitute p\ ot (t") with p\ ot (t) + ff dt"'p H I ot (t"'); however, 
since the terms with derivative of pj ot give contributions 
cither of order w or V, they can be neglected. Thus, we 
can effectively substitute pf^t") with pf^t) in Eq. (15) 
without introducing further approximations. The master 
equation (2) is obtained by rearranging the integration 
limits and the commutators of the last two terms of the 
resulting equation. 



for a systematic accounting of the higher order terms. 
We have two different perturbations, wi(t) and Vi(t), 
and we perform an expansion taking into account all the 
terms up to the order wVV. According to [2] the master 
equation reads 



Pi(t) 



f 



(16) 



Alternative Derivation of the master equation 
(diagrammatic approach) 

A standard iterative derivation of the master equation 
[1] can be cast into a diagrammatic form [2] which allows 



where Xj is the self-energy given by the sum of all irre- 
ducible diagrams. The irreducible diagrams relevant for 
the order wVV are shown in Fig. 3. This produces the 
following non-Markovian master equation: 



Pi{t) = -i[w I (t),p I (t)}--^Tr E J dt'[V I (t),[V I (t'),pj(t')®PB]] 
+ ^ Trs / dt ' I, d *"[vK*), h/(i")JV>(*'),/5/(t')®PB]] 



(17) 



The three terms in the RHS of (17) correspond to the 
diagrams a), b), and c) of Fig. (3) respectively. 

A naive Markovian approximation would be to substi- 
tute pi(t') by pi(t) in both the second and the third terms 
of the RHS of (17). We, however, note, that a Markovian 
approximation consistent with keeping the third term of 
the RHS of (17) requires a substitution of the following 
expression into the second term of the RHS of (17): 

Pi(t') « ~ P i{t) + i j\i!'[ Wl {t"),pi{t")] ■ (18) 



That is, in making the Markovian approximation in the 
second term of the RHS of (17), we have to take into 
account the correction provided by the first term of the 
RHS of (17). We obtain 



a) ^X^Z 

wi(t) 
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FIG. 3: Three diagrams contributing to the self-energy. The wavy line corresponds to the bath correlators (X(t)X(t')) and 
(X(t')X(t)). The horizontal lines are the pieces of the Keldysh contour. The vertices Vi and wi can be put onto the upper (with 
prefactor —i) or the lower (with prefactor i) Keldysh line. This renders the commutators in Eq. (17). In explicit diagrammatic 
notations of Ref. [2] the diagram a) corresponds to two diagrams, the diagram b) to four, and diagram c) to eight. 



1 /•* 

p T (t) = -i[w I (t),p I (t)]-—Tr E J ^'[Vt^JVt^O.P/W^Pb]] 

>/(*),[«;/ (i"), [Vi (t'),pi(t')® Pe}] 

f dt' f dt" fVj(t), [V T (t'), [ Wl (t"), ~ Pl {t") ® p E }] 
Jo Jt' L 



h 2 



(19) 



Finally, upon making the Markovian approximation in ranging the commutators we obtain the following master 
the third and the fourth terms of the RHS of (20), i.e., by equation 
substituting pi(t') by pi(t) and pi(t") by pi(t) and rear- 



1 /•*_.- 
pi(t) = -i[wi(t),pi{t)] - jpTiE J dt'[V I {t),[V I {t'),pi{t)®PE]] 

>j(t),[MO.ft(0],M*)®P* 
I 



+ ^ E L dt 'L dt " 



(20) 
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